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We examine experimental magnetic susceptibility for CaV4 09 by fitting with fitting 

function ax^'^^iT) + c. The function x'^'^^i'T) is a power series of 1/T and the lowest order 
term is fixed as C/T, where C is the Curie constant as determined by the experimental 3- value 
((?=1.96). Fitting parameters are a, c and expansion coefficients except for the first one in 
X™^^(T). We determine a and c as a ~ 0.73 and c ~ for an experimental sample. We interpret 
Q as the volume fraction of CaV409 in the sample and x™^^(2^) ^ tbe susceptibility for the 
pure CaV40g. The result of a 7^ 1 means that the sample includes nonmagnetic components. 
This interpretation consists with the result of a perturbation theory and a neutron scattering 
experiment. 
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§1. Introduction 

A few years ago, Taniguchi et al.& measured the 
magnetic susceptibility of layered material CaV409 and 
found a spin gap of 107 K by analyzing its temperature 
dependence. Since then CaV409 has been investigated 
as an interesting example of two-dimensional spin sys- 
tems with spin gap. There is much theoretical effort 



to undeci 



m Heisenberg 



was not successful to consistently explain various exper- 
iments. In particular it was very difficult to find the 
values of exchange parameters to reproduce the experi- 
mental roagnetic susceptibility and the value measured 
by ESR.N) 

From the lattice structurejHy) the system is expected to 
be described by the two-dimensional Heisenberg model: 



(1) 



H — ^ Jij Si ■ Sj , 

where Si is the spin at vanadium site i. The exchange 
parameter is nonzero if the pair < i,j > corresponds 
to a link indicated in Fig. 1: i. e. Jij is for an edge- 
sharing plaquette link, for an edge-sharing dimer link, 
Jc for an corner-sharing plaquette link and for an 
corner-sharing dimer link. Almost models which have 
been considered are included in this model as special 
cases. 

We summarize important experimental results as fol- 
lows: 

(a) Temperature dependence of magnetic susceptibility 
shows that the material has _a_apin excitation gap 
and its value is about 110 K.Elcil) Measurements, of 
the NMR relaxation time confirms this value.ErES^ 

(b) ESR measurement precisely determined the g-vahie 



(c) 



as 5 = 1.96, which is rather close to 20' The Curie 
constant is then C = 0.003713 emu/g. 
Neutron scattering experiment is performecLto give 
a dispersion relation for spin excitations.E3 The 
result shows a spin gap at momentum (0,0) and 
its value consists approximately with the values of 
other experiments. Also the analysis of the scatter- 
ing intensity directly suggests Jc>Je- 
Theoretically the consistency between the suscepti- 



Fig. 1. Lattice structure for vanadium spins in a layer of 
CaV409. The Heisenberg model examined in text includes 4 
dependent exchange parameters Je (bold solid line), (solid 
line), Jc (bold dashed line) and (dashed line). 
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bility and the g-value is the issue. JBefore experiments 
(b) and (c) appeared, Gelfand et al.B-' estimated the ex- 
change parameters as Je — J'^ — 190 K by assuming 
Jc — J'c — Jel'2-hy using various expansions. In the esti- 
mation, they determined the g-value ds g — 1.77, which 
is different from g = 1.96 in experiment (b). The val- 
ues for the exchange parameters produce the minimum 
at momentum (tt, tt) in the dispersion relatioiL against 
(0, 0) in experiment (c). In a previous article,Q' we es- 
timated exchange parameters as Je(= J^) ~ 610 K and 
Jc{= J'c) — 150 K with the assumption oi g = 2. The 
estimation comes from the determination of the lowest- 
order coefficient in the high temperature expansion of 
the susceptibility. To confirm and refine this result we 
calculated the high-temperature expansion to the third 
order and carried out the fitting to the experimental sus- 
ceptibUity data with keeping the measured g-value {g — 
1.96).E3 However there is no solution satisfying the ob- 
tained set of equations to determine the exchange param- 
eters. Katoh and Imada examined various possibilities 
by taking account of plural atomic levels. However their 
theory includes discrepancy in the compatible explana- 
tion of the susceptibility and the g-value in experiments. 
These theories shed light on some aspects of the suscepti- 
bility of this material, although no one succeeded to give 
a consistent explanation of all the experimental results. 

In this paper, we overcome this difficulty by taking ac- 
count of the possibility that the experimental sample is 
a mixture of CaV409 and other nonmagnetic materials. 
We only assume that the susceptibility can be generally 
expanded in a power series of l/T. From the coefficient 
of l/T in the expansion formula, we estimate the vol- 
ume fraction a of the CaV409 component in the sample. 
We further show that the value of a is consistent with 
that obtained by exchange parameters determined by a 
neutron scattering experiment and the third order per- 
turbation calculation. 

§2. Determination of volume fraction 

To solve the discrepancy we consider that the suscep- 
tibility for the sample is smaller than the true value by a 
constant factor. This is possible if the sample is a mix- 
ture and includes nonmagnetic components. The con- 
stant factor is then the volume fraction of the magnetic 
component. We ha|te two sets of data for susceotibility 
by Taniguchi et al.El' and by Isobe and Ueda.ElP They 
are shown in Fig. 2. Since the latter set is definite for 
the temperature region under 300 K, we compared the 
two sets below 300 K. As a result, we found that the 
latter is precisely 1.17 times lager than the former over 
all temperatures below 300 K. Hence we infer that the 
experimental susceptibility includes a sample dependent 
factor corresponding to the amount of nonmagnetic com- 
ponents. 

In general, the susceptibility is expanded in a power 
series of l/T and the high temperature behavior is de- 
scribed by the expansion formula. This is true for even 
the present case where the sample may be a mixture. 
We make a fitting to the experimental data x*°* (T) by a 
polynomial of sufficiently large order to determine the 
expansion coefficients. To make the expansion coeffi- 



Fig. 2. Experimp«jtal susceptibilities of Taniguchi ct al.liP and of 
Isobe and UedaEiP. The result of 1.17 times the former is also 
shown. 



Table I. Expansion coefficients for susceptibility 



L 


ao 


0.1 


02 


03 


04 


12 


0.00117 


0.72014 


-1.2709 


1.4170 


-1.0914 


13 


0.00188 


0.71095 


-1.2276 


1.3167 


-0.9596 


14 


0.00039 


0.73013 


-1.3176 


1.5224 


-1.2233 


15 


0.00092 


0.72310 


-1.2830 


1.4379 


-1.1060 


16 


-0.00106 


0.74988 


-1.4173 


1.7711 


-1.5742 


17 


-0.00048 


0.74183 


-1.3752 


1.6604 


-1.4074 


18 


0.00001 


0.73513 


-1.3409 


1.5733 


-1.2814 


19 


0.00043 


0.72941 


-1.3123 


1.5028 


-1.1832 


Av 


0.00041 


0.73007 


-1.3181 


1.5252 


-1.2283 



cients dimensionless we take the expansion parameter as 
X = Tq /T with arbitrary constant Tq of the dimension of 
temperature. We choose it as To = 100 K without spoil- 
ing the generality. The polynomial for fitting function is 
accordingly written as 

Here the expansion coefficients a,„'s are fitting parame- 
ters. The order L of the polynomial is chosen to be suf- 
ficiently large so that a„i's do not change largely when 
L changes by a few integers. 

We carried out the fitting for each L less than 19 and 
obtained the values of expansion coefficients. Results for 
L = 12 to 19 are shown in Table I. Averages over L 
for 12 to 19 are also shown in the last low. Viewing 
this table we see that the values for gq and ai have only 
small fluctuation and are relatively reliable. In contrast 
the values for seem to largely fluctuate and include 
an amount of fltting error. The fltted curve is shown in 
Fig. 3 together with x*'°HT)- 



Magnetic Susceptibility for CaV409 



3 



Fig. 3. The fitting function for the optimal fitting to the exper- 
imerijtial susceptibility. Shadowed circles are experimental data 
X*°* which seems to form a line by overlap in a high temper- 
ature region. 



The polynomial ^ is rewritten as 
C 

/(x) = — [ao + aia;(l + a^jX + CgX^ H )] (3) 

-to 

with = Qi/ai for i > 2. Correspondingly we analyze 
the experimental susceptibility x*°'(^) the following 
form: 

x*°*(r) = 7^ + «x"^^^(r), (4) 

Jo 

where j — and a — ai. Then x"^^^(2^) is approxi- 
mately represented as 

^^--^{T)^^{l + A,x + A,x' + ---) (5) 

with Ai = a'i^i. Hence we estimate the values of 7 and 
a as 

7 ~ 0.00, a ~ 0.73. (6) 

For an ideal localized spin system, no constant term 
appears in the magnetic susceptibility. The constant 
term of 7 in eq. comes from other than the electronic 
spins if it exists. Candidates for the contribution are 
the Van Vleck term and the diamagnetic term. However 
these terms have been already subtracted in the present 
data of eq. The present result, 7 = 0, confirms 

that the original subtraction is precise. We note that 
the present method uses only the general property that 
the susceptibility is expanded in power series of 1/T. In 
contrast the method for the original subtraction explic- 
itly uses the facts that the system is a spin system and 
is two-dimensional. 

If the sample consists only of CaV40g, the observed 
susceptibility must be x^^'^iT) instead of x*°*(r). This 
is because the Curie constant C is fixed as C = 0.003713 
emu/g due to experiment (c) and the constant a must 



be unity. The result of a < 1 means that the sample 
includes nonmagnetic materials, which do not have spin 
degree of freedom. Precisely the sample is a mixture 
consisting of 73 % magnetic and 27 % nonmagnetic ma- 
terials. Thus a is the volume fraction of the magnetic 
part in the whole sample. We infer that the nonmagnetic 
part is due to incomplete chemical reactions and hence 
CaO and VO2 are candidates. We believe that this is 
the reason why the effective g-value in some past the- 
ories are much smaller than the observed value by the 
ESR measurement. 

§3. Consistency with the Perturbation Theory 

As mentioned in experiment (c), Kodama et al. ob- 
tained the dispersion relation of, a spin excitation for 
CaV409 by neutron scattering^' They also theoreti- 
cally calculated a dispersion formula by a perturbation of 
the second order for the Hamiltonian (|l|). By fitting the 
formula to the experimental dispersion, they determined 
the exchange parameters. The resultant values approxi- 
mately reproduced the experim(3ii|tal dispersion relation. 
Recently, Fukumoto and Oguchitj' recalculated and cor- 
rected the second order perturbation formula of the dis- 
persion relation. They further extended the perturbation 
formula up to the third order. In the third order pertur- 
bation, their result is given as Jc = 162 K, Jc = 20 K 
and J'g = Je = 79 K. Using the third-order values for the 
exchange parameters, we numerically diagonalized the 
Hamiltonian and obtained energy levels. The diago- 
nalization was done for lattices with N — 8 and 16 (A'': 
the number of lattice sites) under the periodic boundary 
condition. We then calculated the magnetic susceptibil- 
ity x""™(2^) by the energy levels. The results are shown 
in Fig. 4. The curve for 8 sites is close to that for 16 
sites. We extrapolate the results to the infinite system 
size by assuming the system size dependence 1/iV. We 
plotted the values of 0.7 times the extrapolated data in 
Fig. 4. The resultant curve is close to the experimental 
susceptibility and shows that the volume fraction is given 
as a' 0.7. This result consists with the interpretation 
obtained from the analysis of the magnetic susceptibility. 

The difference between a and a' is not serious since 
some errors may be included in the processes determining 
them. There is a possible error for a in the fittingj-of 
a polynomial (||) to the experimental data x'°*(T).Ej' 
Further there is a possible error in a' accompanied with 
the truncation of the perturbation series to the third 
order. 

§4. Summary 

We examined the experimentally obtained magnetic 
susceptibility x'°*(^) for CaV409. By using the fact 
that the susceptibility is generally expanded in the power 
series of l/T, we showed that the sample includes non- 
magnetic components. We estimated the volume frac- 
tion of CaV409 as a ~ 0.73. The existence of nonmag- 
netic components enable us to consistently understand 
all the experimental results: the magnetic susceptibil- 
ity, the g-vahie by ESR and the dispersion relation by 
neutron scattering. 
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[22] T. Ohama, H. Yasuoka, M. Isobe and Y. Ueda: J. Phys. Soc. 
Jpn. 66 (1997) 23. 

[23] To estimate the accuracy of the fitting method in §4, we ap- 
pUed the same fitting to the data of the numerical diagonal- 
ization {N = 16) instead of the real experimental data. We 
found that approximately 2 % error is involved in the esti- 
mated coefficient for the first order term, corresponding to 
the volume fraction. 



Fig. 4. Susceptibility obtained by the numerical diagonalization 
of the Hamiltonian (^. Dash-dotted (dashed) line is for N = 8 
(N = 16) lattice. The solid line is for 0.7 times the large N limit 
of these numerical results. Experimental data are also shown. 
The values of the exchange r^rameters in the Hamiltonian are 
what Fukumoto and Oguchit£P determined. 
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